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We present a calculation of protostellar disk formation and evolution in which 
gaseous clumps (essentially, the first Larson cores formed via disk fragmenta- 
tion) are ejected from the disk during the early stage of evolution. This is a 
universal process related to the phenomenon of ejection in multiple systems of 



CLi! point masses. However, it occurs in our model entirely due to the interaction of 

O ! compact, gravitationally-bound gaseous clumps and is free from the smoothing- 

^ I length uncertainty that is characteristic of models using sink particles. Clumps 

c^, | that survive ejection span a mass range of 0.08-0.35 Mq, and have ejection ve- 

locities 0.8 ±0.35 km s~^, which are several times greater than the escape speed. 
^ I We suggest that, upon contraction, these clumps can form substellar or low- mass 

!^ I stellar objects with notable disks, or even close-separation very-low-mass bina- 

CS I ries. In this hybrid scenario, allowing for ejection of clumps rather than finished 

protostars/proto-brown-dwarfs, disk formation and the low velocity dispersion 



Q I of low-mass objects are naturally explained, while it is also consistent with the 

CN I observation of isolated low-mass clumps that are ejection products. We conclude 

that clump ejection and the formation of isolated low mass stellar and substellar 
objects is a common occurrence, with important implications for understanding 
the initial mass function, the brown dwarf desert, and the formation of stars in 
all environments and epochs. 
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mass, brown dwarfs 
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Introduction 



The formation of brown dwarfs (BDs) is usually attributed to either of two distinct 
mechanisms. One model is the direct collapse of a very low mass cloud fragment whose 
mass straddles the substellar mass limit. Formation mechanisms for such low mass (and 
consequently high density) cores include colliding flows in a turbulent magnetic medium 
( jPadoan &: Nordlundl 120041 ) and photoerosion of more massive cores by the ionizing UV 
radiation from nearby OB stars (jWhitworth &: Zinneckerl |2004| ). Since the collapse of a 
prestellar core has been shown by numerous studies to lead to a very focused collapse with 
a power-law density profile, the probability of fragmentation during the runaway collapse 
phase (i.e., before the formation of a central hydrostatic object) is low. This leads to the 
scenario that the stellar and substellar mass function is simply a scaled-down version of the 
core mass function. 

However, during the protostellar phase (i.e., after the formation of a central hydro- 
static object), and unlike in the prestellar phase, the relative scaling of flow variables in 
the p resence of angular moment um can lead to the formation of a centrifugally-balanced 
disk (JBasu fc Mouschoviad 119951 ). Such a disk has a high density and evolves on a time 
scale much longer than its dynamical time. Therefore, the protostellar disk has time for 
fragmentation if conditions are suitable, e.g., the criter ia for gravitational instability and 
fragmentation are satisfied ( lToomrelll981t lGammiell200ll ). Furthermore, fragmentation may 
be induced by clos e encounters with members of a stellar cluster in otherwise gravitationally 
stable disks (e.g., iThies et al.ll2010[ ). This leads to the second model: that disk fragmen- 



tation is responsible for the formation of low mass and substellar objects in the vicinity of 
a forming star, even if the mass of the parent core is far greater than the substellar limit. 
In this scenario, free-floating BDs would be the result of ejections from the disk, which is a 
plausible outcome based on extrapolation from three-body calculations of interacting point 
masse s and can account for pe c uliar proper ties of BDs as compared to those of low-mass 
stars (JKroupa fc Bouvierl 120031 ). iBatd (120091 ) obtained the BD ejection from SPH modeling 
of molecular cloud fragmentation resolved from pc scales down to the opacity limit that 
occurs on AU scales. In that model, multiple collapsing fragments each form a disk. Some 
disks undergo fragmentation, and low-mass objects, treated as point masses in t he SPH ap- 
proximation, can be ejected from multiple systems. In the models of lBatd ( 120091 ). a sink cell 
is introduced on AU sc ales and gas is cleared away frorn a region around them of size up 
to 10 AU. Furthermore, IStamatellos fc WhitworthI ( l2009l ) have modeled the formation and 
ejection of BDs after disk fragmentation, also using SPH and sink cells, and starting from an 
initial condition of a pre-existing massive disk that has a mass equal to that of the central 
star. 
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In this paper, we pursue models that solve the self-consistent formation of a centrifugal 
disk from the collapse of a larger cloud core, employing the thin-disk approximation. Compu- 
tational efficiency allows us to run a large number of models, with various initial conditions 
of cloud mass and angular momentum, as well as study the long-term evolution of disks. We 
employ a central sink cell of size 6 AU, representing a central star and some circumstellar 
material, but the surrounding disk that forms is followed without recourse to sinks, even 
when fragmentation occurs. Our main result is that ejection of protostellar/protobrown 
dwarf mass clumps occurs through the interaction of multiple gravitationally-bound gaseous 
clumps, which are essentially the first Larson cores formed via gravitational fragmentation 
of protostellar disks. 

Our results support a new hybrid paradigm of BD and low mass star formation by 
disk fragmentation followed by clump ejection (rather than by ejection of finished low-mass 
stars/brown dwarfs). This scenario naturally accounts for the presence of disks around BDs 
and the relatively low ejection speeds in this scenario can account for the observed velocity 
dispersion and physical location of BDs relative to YSOs. It is also consistent with the so- 
called brown dwarf dese rt, i.e., the lack of BD companions to primary star s at distances less 
than several tens of AU ( McCarthy fc Zuckermanll2004l : iMarcy et al. 20051). A disk fragmen- 
tatio n origin for BDs can explain the bro wn dwarf dese rt (see, e.g.. lStamatellos fc Whitworth 
201ll ) since cooling constraints in disks ( iRafikovl 120071 ) favor fragmentation at radii > 50 AU 
from the primary object. Furthermore, fragments that do form at closer range or are torqued 
inward from wider orbi ts are often subject to a run away inward migration toward merger 
with the central object (jVorobyov fc Basull2005l . |2006[ ). so that low mass gaseous companions 
are not typically expected to settle into stable inner orbits. However, it has been suggested 
that some of the eventually tidally disrupted clumps may have enough time to build up 
solid cores. These may then be released into inner orb its after the disruption, proyiding 



thermally-processed sol i ds or even rocky planetary cores ( jBoley et al.ll2010t iNayakshinl 12010 



Cha &: Nayakshinll2011t IVorobyovl 1201 ll ) . This scenario may account for planet formation in 
the inner disk but presumably not something as massive as a BD. Detailed simulations of 
this possible mechanism remain to be worked out. 

In contrast, the direct core collapse paradigm for BD formation also suffers from the 
need to explain very low mass gravitationally bound objects that are one to two orders 
of magnitude less than the ambient Jeans mass in molecular clouds. This can be ac- 
complished with very high bounding pressure, wheth er due to very strong magnetic fields 



(iBasu fc Mouschoviaslll995l ) or high turbulent pressure (jPadoan &: Nordlundll2002 



Chabrier fc Hennebelle 



201ll ). Both the ideas of extremely high magnetic or turbulent bounding pressu re on a frag 



ment suffer from a lack of direct observational evidence at the moment (see Andre et al. 



20091 ). Turbulent pressure may also not be present in such a focused manner as to act like 
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an isotropic bounding pressure. 



Model 



Our model and method of solution are the ones pre sented in detail bvlyorobvov fc Basu 



(120 lOl ). which is an update to the model presented by IVorobyov fc Basul (120061 ) . We solve 
the mass, momentum, and energy transport equations in the thin-disk approximation, us- 
ing vertically-integrated quantities. Disk self-gravity is calculated using a two-dimensional 
Fourier convolution theorem for polar coordinates. The energy equation includes compres- 
sional heating, heating due to stellar and background irradiation, radiative cooling from 
the disk surfaces, heating due to turbulent viscosity (parameterized via the usual Shakura 
& Sunyaev prescription with a = 0.005), and a diffusion approximation to link the effec- 
tive surface temperature of the disk with the midplane t e mpera ture. Frequency-integrated 
opacities are adopted from the calculations of JBell fc LinI (119941 ) and a smooth transition is 
introduced between the optically thin and optically thick regimes. In the opaque regime, 
gravitationally-bound clumps that form in the disk are essentially " f irst hydrostatic cores" 
or "fi rst Larson cores" in the parlance of star formation (JLarsonlll969l : iMasunaga fc Inutsuka 
2000l ). A "second collapse" down to stellar dimensions can take place when the central 
temperature reaches ~ 2000 K and H2 is dissociated. This possibility is not allowed in our 
energy equation so that we may continue to track the evolution without recourse to sink 
cellq^, however clumps in the disk typically do not reach this temperature before disrup- 
tion or ejection (see Section 4). The stellar irradiation is based on a luminosity that is a 
combination of accretio n luminosity and photosphe ric luminosity calculated from the pre- 
main-sequence tracks of lD'Antona fc Mazitellil (jl997| l. The pressure is related to the internal 
energy through an ideal gas equation of state. 

Models presented in this paper are run on a polar coordinate (r, (p) grid with 512 x 512 
zones. A central sink cell of radius 6 AU is employed and the potential of a central point 
mass is added to the disk self-gravity once a central hydrostatic core is formed. The radial 
points are logarithmically spaced, with the innermost cell outside the central sink having size 
0.07-0.1 AU depending on the cloud core size (i.e., the radius of the computational region). 
The latter varies in the 0.025-0.12 pc (5000-24000 AU) limits. The radial and azimuthal 
resolution are about 1 AU at a radius of 100 AU. The inner and outer boundary conditions 



"'^ Resolving second cores with a typical size of several stellar radii, and thus avoiding the use of sink 
particles, is a formidable task for grid-based codes designed to study the global evolution of protostellar disks 
as it requires a local spatial resolution < 0.01 AU. 
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are set to allow for free outflow from the computational domain. The initial conditions for 
the prestellar core correspond to profiles of column density E and angular velocity Q of the 
form 



\/r'^ + r\ 



Q 



21]n 



r ) 




(2) 



These profiles have a small near-uniform central region of size r^ and then transition to an 
r~^ profile; they are representative of a wide class of observations and theoretical models 
flBasulll997l : lAndre et al.ll2009l : iDapp fc Basull2009l ). Our previous papers have also shown 
that the qualitative features of the ultimate centrifugal disk evolution are independent of the 
specific profiles adopted for the initial prestellar core. Our reference model presented in this 
paper has rotational support in the initial state characterized by the dimensionless parameter 
1] = ^r\j(?^ = 1.4 X 10~^, where Cg = 0.188 km s^^ (corresponding to a temperature 
T = 10 K) in the initial state but varies spatially and temporally during the evolution. 
Alternately, the ratio of rotational energy to the magnitude of gravitational potential energy 
is /3 = 1.3 X 10"^. Table 1 shows a list of parameters and some results from models we have 
run for this study. The reference model of our study (model 5) has tq = 2006 AU, a total 
core mass M^ = 0.92 Mq, and /3 = 1.3 x 10"^ 



3. Results 



Our recent numerical hydrodynamics simulations with realistic cooling have shown that 
fragment formation is enhanced at radii > 50 AU (in agreement with many previous studies 
on disk fragmentation), and that e pisodic accretion due to clump formation and inspiralling 
of the clumps is a robust result fIVorobyov fc Basull2010l ). In this paper, we report on a 
qualitatively new kind of result — an occasional ejection of fragments from the disk into the 
intracluster medium. Long integration times (~ 1 Myr) need to be run in order to capture 
the occasional ejection events. Below, we discuss in detail the results for our prototype 
model. 

Figure 1 shows a sequence of column density images for our reference model (model 5), 
at various times after the formation of a central object, and in a region of size 2400 AU 
on each side. The circumstellar disk forms at t ~ 0.01 Myr and already at t = 0.05 Myr 
it undergoes fragmentation. Within the first ^ 0.25 Myr, multiple fragments are formed 
in the relatively massive disk, at distances of > 50 AU and < few hundred AU. They 
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Table 1: Models showing ejection and properties of ejected clumps 



Model 


/3 


^0 


'^o 


Me 


Ms 


Meject 


''eject 


"Reject 


"^esc 


It 


0.23 X 10-2 


0.48 


3940 


1.78 


1.03 


0.13 


0.52 


0.38 


0.34 


2* 


0.56 X 10"2 


1.14 


2400 


1.08 


0.68 


0.15 


0.34 


0.77 


0.31 


3t 


0.56 X 10-2 


0.72 


3770 


1.69 


0.94 


0.06 


0.43 


0.6^ 


0.33 


4t 


0.56 X 10-2 


0.88 


3085 


1.38 


0.82 


0.05 


0.70 


0.5* 


0.36 


5 (ref.) 


1.3 X 10-2 


2.0 


2006 


0.92 


0.47 


0.15 


0.30 


0.91 


0.31 


6 


1.3 X 10-2 


1.2 


3430 


1.54 


0.45/0.69 


0.2/0.2 


0.24/0.65 


0.88/1.20 


0.26/0.27 


7 


1.3 X 10-2 


1.0 


4115 


1.85 


0.86/0.87 


0.08/0.2 


0.73/0.83 


0.76/0.52 


0.42/0.29 


7t 


- 


- 


- 


- 


0.73 


0.11 


0.45 


0.7* 


0.45 


8 


1.3 X 10-2 


3.4 


1200 


0.54 


0.31 


0.1 


0.18 


1.0 


0.32 


9t 


1.3 X 10-2 


1.5 


2780 


1.25 


0.49 


0.2 


0.22 


0.45* 


0.29 


io*t 


2.24 X 10-2 


2.9 


1885 


0.85 


0.46 


0.08 


0.54 


0.44 


0.28 


ir 


2.24 X 10-2 


1.4 


3940 


1.77 


0.51 


0.34 


0.33 


0.45 


0.32 


lit 


- 


- 


- 


- 


0.56 


0.09 


0.45 


0.62* 


0.38 























Note. — All distances are in AU, velocities in km s ^, angular velocities in km s ^ pc ^, time in Myr, 
and masses in Mq. 
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are generally torqued inward through gravita tional interaction with trailing spiral arms, as 



first found by lVorobyov fc Basul (120051 . l2006l ). Others located at large radii may eventually 
disperse. However, under sometime favorable conditions, a clump within a multi-clump 
environment can be ejected through many-body interaction. The ejection is also aided by 
the non-axisymmetric potential of the relatively massive disk. The arrow in the image in the 
second row of the middle column points to a clump that is subsequently ejected from the 
system. The velocity of the ejected clump (at the moment when it leaves the computational 
boundary at r = 12000 AU) is about 0.9 km s~^, which is a factor of 3 greater than the 
escape velocity tiosc = (2GMsd/'r)^'^, where the total mass of the central star and its disk is 
Msd = 0.6 AIq. This means that the clump will truly be lost to the system. The ejection 
event is transient and the clump is not seen in the subsequent image at t = 0.27 Myr. The 
total mass of the clump, calculated as the mass passing the computational boundary during 
the ejection event, is 0.15 Mq. It should be noted, however, that this value includes not 
only the compact core but also a diffuse envelope and even fragments of spiral arms (see 
Fig. 2 below). We speculate that upon contraction this clump may form a substellar object, 
given that a significant fraction of mass remains in the disk until it is ejected due to outflows 
and/or dispersed due to photoevaporation. No further major fragmentation episodes in the 
protostellar disk are seen after t = 0.27 Myr, likely due to a sharp drop in the total disk 
mass caused by the ejection. The disk then gradually evolves toward a nearly axisymmetric 
state. 

To illustrate the ejection mechanism in more detail. Figure 2 shows a sequence of column 
density images on a much larger spatial scale of the simulation box, 20000 AU on each side, 
and on a much shorter time scale, t =0.26-0.31 Myr. The arrows in the images identify the 
same clump that was identified in Figure 1, and show its ejection from the disk and indeed 
the entire simulation box. As the clump moves nearly radially outward, it grows in size and 
diffuses somewhat but this is an artefact of the decreasing resolution of our numerical grid 
at larger radii. Future higher resolution simulations will be needed to verify the evolution of 
the clumps after they undergo an ejection event. 

Figure 3 shows the mass accretion rate onto the central sink as a function of time since 
the beginning of numerical simulations. The time t = 0.12 Myr represents the moment that 
a central object is formed. After a brief period of smooth accretion, a disk forms outside the 
sink cell, and the subsequent accretion is highly variable, with the highest rates exceeding 
10~^ Mq yr~^ and the lowest rates plunging below 10~^ Mq yr~^. Interestingly, there is a 
late burst at 0.38 Myr (0.26 Myr after the formation of the central star) that is greater in 
magnitude than all previous bursts. This burst corresponds to the inward flow of material 
that accompanies the ejection of the clump. The effects are twinned due to the need for 
overall angular momentum conservation — one clump is ejected while the other one is driven 



0.5 




-500 500 

Radial distance (kpc) 



Fig. 1. — Gas surface density distribution (g cm"^, log units) in the reference model at 
several time instances after the formation of a central star. The box size is 2400 AU on 
each side, and represents a small subregion of the overall computational domain. An arrow 
identifies a clump that is ejected from the inner region soon afterward. 
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Fig. 2. — Gas surface density distribution (g cm"^, log units) in the reference model at 
several time instances after the formation of a central star. The box size is 20000 AU on 
each side, and represents nearly the full extent of the computational domain. Arrows identify 
a clump that is ejected from the system after a multi-body interaction within the centrifugal 
disk. 
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into the center. After this last major burst, the disk settles down to a more quiescent phase of 
evolution with a much lower accretion variability, characterized by flickering in the accretion 
rate but no bursts, and a gradual decline in accretion rate. The residual accretion rate of a 
few X 10~^ Mp; yr~-^ at 1.0 Myr is consistent with observed accretion rates for T Tauri stars 



[e.g. 



Calvet et al.ll2004l ). This behavior agrees well with a gradual gravitational stabilization 
of the disk that is caused the ejection event and the associated substantial mass loss. From 
this time onward, viscous torques start to dominate gravitational torques in the disk mass 



transport (IVorobyov &: Basull2009l ). 



Figure 4 illustrates the evolution of mass in various components of the system, as well 
as in the overall simulation region. The total mass in the system suffers a drop as the ejected 
clump exits through the outflow outer boundary applied to the simulation region. The mass 
loss equals the clump mass: ^0.15 Mq. The clump maintains its form all the way to the edge 
of our simulation box notwithstanding a constantly decreasing resolution on the logarithmic 
grid, and is gravitationally bound, with the magnitude of its gravitational potential energy 
being greater than the sum of its rotational and thermal energies. Hence, we expect it to 
form a low mass star or substellar object depending on the efficiency with which it converts 
into a collapsed object. Figure 4 also shows a notable sharp drop in the disk mass and an 
equivalent increase in the envelope mass at t = 0.26 Myr. These effects are the manifestation 
of the ejection event when the clump leaves the parent disk and joins the envelope. At about 
the same time, the stellar mass increases sharply due to consumption of the twin clump that 
has lost its angular momentum due to the many-body interaction. 

We have considered 20 models that exhibit disk fragmentation. The corresponding 
initial prestellar core masses lie in the range Mc = 0.3 — 1.8 M©, and ratios of rotational to 
gravitational energy are in the range /3 = (0.25 — 2.3) x 10"^. Only 11 models have shown 
ejection events and the parameters of these models are described in Table [TJ A star symbol 
next to the model number means that the ejected clump was initially a binary clump, and 
a dagger symbol means that the clump is smeared out before reaching the computational 
boundary. Both these phenomena are discussed later in this section. The last five columns 
present the mass of the central star at the moment of ejection Mg, the mass Mgject of the 
ejected clump, the time teject after the formation of the star when the ejection event takes 
place, the ejection velocity Reject, and the escape velocity Wesc, respectively. Both Weject and 
Wesc are calculated at the time when the clump passes through the outer computational 
boundary. A diamond symbol indicates that the corresponding value is calculated at a 
smaller distance, since the ejected clump has been smeared out before reaching the outer 
computational boundary. When more than one ejection event is present in the same model, 
the different ejection speeds are separated by a slash. In all cases, Reject > ^esc, indicating 
that the ejected clumps will be lost to the system and become freely-floating objects. The 
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mean and median ejection velocities in our models are 0.7 km s~^ and 0.6 km s~^, and the 
minimum and maximum values are 0.38 km s~^ and 1.2 km s~^, respectively. These are 
consistent with the velocities of low mass stars and BDs in the Chamaeleon I and Taurus 



star- forming regions (e. g. . I Joergensll2006l ) . There is practically no dependence of the ejection 



velocities on the mass of the ejected clumps. 

Figure [5] presents the total mass in the system and illustrates the ejection events in 
11 models of Tabled] A sharp decrease in the total mass indicates the time when the clump 
escapes the computational region. However, approximately half of the models (marked 
with a dagger symbol in Table [1]) demonstrate a rather smooth decline in the total mass, 
extended over (1 — 2) x 10^ yr. In these models, clumps have dispersed before reaching the 
outer computational boundary. They may become gravitationally unbound due to strong 
tidal torques during the close encounters that lead to ejection, and/or due to decreasing 
numerical resolution on a logarithmically spaced grid as they propagate radially outward. 

Table [1] demonstrates that some of the ejected clump masses Meject are below the sub- 
stellar mass limit, giving us more confidence to propose this as a potential BD formation 
mechanism. Note that because this process occurs due to the interactions of gaseous clumps 
and a gaseous disk, it is somewhat different from the scenario of ejections by multi-body 
interactions of point masses, i.e., objects that have already fully collapsed down to compact 
equilibrium configurations. 

Analysis of the ejection times in Table [T] and Fig. E] shows that most ejection events 
take place between 0.2 Myr and 0.8 Myr after the formation of the central star and no 
ejections are seen after t = 1 Myr. This is due to the fact that gravitational instability and 
fragmentation are limited exclusively to the early stages of disk evolution. 

It is interesting to note that models 6, 7, and 11 experience multiple ejection events. 
The same may may be true for model 10 but the ejected masses at t ^ 0.25 Myr and t ~ 0.4 
Myr are too small for the events to be resolved. In addition, in models 2, 10, and 11 (marked 
with the star symbol), the ejected clumps were initially close binaries with separations of 
order a few AU, but later merged as they propagated toward the outer boundary, most likely 
due to decreasing resolution at large radii. The possibility of ejection of close binary clumps 
is a particularly appealing feature of the clump ejection scena rio for BD formation si nce the 



separation of binary BDs seems to peak at around 1 — 4 AU (jWhitworth et al.l 120071 ). 



Fig. E] shows the P — M^ phase space covered in our modeling. The dark-shaded area 
highlights the region where both disk fragmentation and clump ejection are observed, while 
the light-shaded area defines the region where only disk fragmentation takes place. There 
is a increasing tendency for ejections to take place as /3 and Mc increase. Models that fall 
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into the light-shaded domain usually form just one fragment at a time, in which case the 
ejection mechanism cannot work even in principle. Furthermore, not all models that are 
in the ejection domain may undergo ejection events due to the highly stochastic nature of 
this mechanism. This means that the dark-shaded area outlines the region where ejection 
is possible in principle, but in reality it might not be observed if a chance arrangement of 
clumps would not favour their close encounters. The robustness of Figure [6] is also affected 
by a rather coarse grid of models explored in the paper. For instance, we have considered 
only four values of f3 and about ten models for each /3. The uncertainty associated with a 
small number of models is illustrated in Figure [6]by arrows showing where the fragmentation 
(and ejection) regions could potentially extend, had we considered a finer grid of models. 
We restricted our modeling with Mc < 2.0 Mq and ratios (3 < 0.025 due to the fact that 
models with higher initial core masses and angular momenta tend to form wide-separation 
(> 10 AU) binary/multiple stellar systems for which our numerical code is not well suited. 
Disks around such systems may be strongly truncated due to tidal effects, which may impede 
disk fragmentation and ejection. On the other hand, close encounte rs with other mern bers 
of a stellar cluster may promote disk fragmentation and ejection (JThies et al.l l2010l ). A 
dedicated study using suitable numerical codes is needed to address these issues. We also 
note that a steep cutoff of the fragmentation region at small /3 < 0.001 is a physical effect 
related to the fact that such models form di sks of too low a mass and size for gravitational 
fragmentation to take place fjVorobyovl 1201 ll ) . 



4. Discussion and Conclusions 



We find that high resolution models with a realistic treatment of heating and cooling 
and in which a disk forms consistently from prestellar collapse, can sometimes lead to the 
ejection of gaseous clumps from the disk and eventually from the system (defined by the 
extent of the original prestellar core). These are the first confirmed examples of ejections 
due to interactions of compact gaseous clumps and their parent gas disks and host star, as 
opposed to being caused b y the i nterac tion of point particles. Previous numer ical simulations 
of ejections from disks by iBatd ( l2009l ) and IStamatellos fc Whitworthl (120091 ) have used sink 
particles with a smoothing length as a free parameter. While the smoothing length can 
in principle be used to mimic the size of a clump, our numerical simulations are free from 
such uncertainty and the size of the clump is instead determined self-consistently by the 
disk thermodynamics. We show explicitly that the clumps formed via disk fragmentation 
need not contract to stellar densities in order to be ejected into the intracluster medium. 
Moreover, some of the clumps may be tidally disrupted during the ejection event and disperse 
as they leave the system, a phenomenon that is impossible to capture using sink particles. 
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Fig. 3. — Mass accretion rate onto the star as a function of the time elapsed since the 
beginning of the collapse in the reference model. There is a particularly large spike in mass 
accretion rate coinciding with the ejection of the clump at 0.38 Myr (0.26 Myr after the 
formation of the central object). 
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Fig. 4. — Time evolution of the masses of the disk, star, envelope (including ejected com- 
panions when present) and total mass as a function of time since the formation of the star. 
Note the large drop in disk mass and concurrent growth of envelope mass at about 0.26 Myr, 
the time that the clump is ejected from the disk. 
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In fact, most of the clumps in the disk never achieve tempe ratures required for a second 



collapse before they a re disrupted, as pointed out previously (lNayakshinll2010l : iBolev et al. 



2OIOI : I Vorobyovl |20 1 ll ) . For instance, figure 2 (bottom- left panel) in IVorobyovl (120 ll[ ) shows 
that most of the clumps have maximum temperatures below 1000 K and only one clump 
achieves a temperature higher than 2000 K before being disrupted or torqued onto the star. 
This makes the phase of interaction of multiple first cores, as studied in this paper, all the 
more relevant to disk evolution. 

An important argument usually in favor of t he direct collapse sce nario is the observation 
of some very low mass isolated proto-BD clumps fJLuhman et al.ll2007l ). However, by allowing 
for the possibility of ejected clumps rather than finished BDs, our hybrid scenario also allows 
for this possibility. Clumps that survive can cool and contract to form a (sub-)stellar core 
and a disk, thus also explaining the presence of accretion disks around BDs. Furthermore, 
the ejection speeds of the clumps are < 1 km s^^ ?« 1 pc Myr~^, reflecting the fact that 
ejections take place in a much more diffuse situation than interacting point masses in close 
proximity, i.e., with impact parameters of tens of AU rather than several radii of finished BDs 
or low mass stars. With these ejection speeds, the velocity dispersion of BDs is not expected 
to differ much from that of YSOs and they should remain spatially co-located within cluster 
forming clumps of pc scale for at least the typical 1-2 Myr age of YSO clusters. Our models 
can also explain qualitatively why the BDs are less numerous than stars, since the ejections 
only occur for some region of parameter space. While a quantitative prediction of BD fraction 
remains difficult without more insight into core initial conditions, we can use a stellar IMF 
to obtain some estimates. Assuming that about half the mass of each clump accretes onto 
a star (i.e., 50% formation efficiency), that only clumps with m ass M^ > 0.8 Mq undergo 
ejection events, and using the Kroupa stellar IMF (JKroupall200ll ) with lower and upper cutoff 
masses 0.08 Mq and 100 Mq, respectively, we estimate that the number of ejected clumps is 
about one for every 10 stars. Here we take into account that clumps survive ejection in only 
about half the realizations, but this decrease is mostly offset by models showing multiple 
ejection events. The net ejection rate is then about one clump per model. T his estimate is 



consi stent with the estimated empirical ratio of 5-10 between stars and BDs (JLuhman et al. 



20071 ). We must keep in mind that our estimate is rather tentative and we expect it to 
increase somewhat if more models are run at higher resolution. There are also hints in our 
simulations that ejected clumps with substantial rotation may form close binary objects, in 
agreement with the prevalence of close binaries in BD-BD systems. All in all, our hybrid 
scenario agrees, at l east qualitatively, with nearly all of the key observational features (see 
Luhman et al.ll2007l ) of BDs. The only observation not easily accounted for are the wide BD 
binaries, but these are very rare and they may indeed form from direct collapse since the 
total mass required is likely not substellar. 
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Our results reveal a fundamental process that is intrinsic to gaseous disks with moderate 
to high levels of angular momentum. This clump ejection mechanism should be robust to 
environmental variations of metallicity and temperature, as it depends more directly on 
angular momentum and multi-body dynamics. Therefore, we believe that low mass stars 
and substellar objects should form in a variety of situations in different cosmic locales and 
epochs. 
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